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DEUTERIUM  DIFFUSION  IN  A  SOIL-WATER-ICE  MIXTURE 

by 

Joseph  L.  Oliphant  and  Allen  R.  Tice 

INTRODUCTION 

It  has  been  suggested  that  deuterated  water  could  be  used  as  a  tracer 
In  experiments  aiming  to  clarify  the  mechanisms  of  water  and  Ice  migration 
in  frozen  soils.  Arnason  (1969)  has  shown  that  for  summer  water  percola¬ 
tion  In  a  temperate  glacier  in  Iceland,  equilibrium  of  deuterium  concentra¬ 
tions  between  the  ice  and  water  Is  closely  approximated.  Arnason  also 
showed  in  a  3-hour  laboratory  experiment  with  melting  snow  that  equilibrium 
was  approached  but  not  reached.  To  use  deuterium  oxide  as  a  tracer  for  wa¬ 
ter  migration  in  partially  frozen  soil,  the  processes  and  rates  of  equili¬ 
bration  of  deuterium  between  frozen  and  unfrozen  phases  need  to  be  under¬ 
stood  and  quantitatively  known. 

The  processes  by  which  equilibration  between  phases  containing  deuter¬ 
ium  take  place  include  diffusive  and  convective  movement  of  water  contain¬ 
ing  deuterium  and  protium  in  the  liquid  phase  and  diffusive  movement  in  the 
solid  phase.  There  also  may  be  some  resistance  to  movement  across  the 
phase  boundaries  that  needs  to  be  considered.  The  diffusion  coefficient 
for  deuterated  water  in  the  liquid  phase  is  on  the  order  of  10“  cm  /s 
(Franks  1982)  and  that  for  deuterium  in  ice  is  about  10“  cm  /s  (Franks 
1972).  The  relatively  low  but  finite  diffusion  rate  in  ice  can  cause,  un¬ 
der  some  conditions,  equilibrium  between  the  ice  and  water  but  under  other 
conditions  this  equilibrium  will  not  be  attained.  The  following  experiment 
was  designed  to  test  the  rate  of  equilibration  of  deuterated  water  between 
the  ice  and  liquid  water  phases  in  frozen  soil. 

EXPERIMENTAL 

Oven  dry  raorin  clay  was  mixed  with  about  IX  by  weight  deuterium  oxide 
(DjO)  and  then  compacted  in  stainless  steel  tubes  of  15-mm  l.d.  to  a  thick¬ 
ness  of  about  5  mm.  About  20  specimens  were  prepared,  placed  in  a  cold  en¬ 
vironment  (-20°C)  and  frozen.  Following  freezing,  each  specimen  or  wafer 


was  sanded  to  a  uniform  thickness  of  2.5  ±  0.1  mm  with  Carborundum  cloth 


and  fine  sandpaper.  Each  sanded  wafer  was  then  stored  in  a  small  airtight 
glass  bottle  to  minimize  sublimation. 

Cylindrical  single  crystals  of  ice  of  15-mm  o.d.  and  about  AO  mm  long 
were  prepared  from  deionized,  distilled  water  by  use  of  a  dye  extrusion 
technique.  These  were  sawed  into  5-mm  thick  sections,  then  sanded  to  a 
thickness  of  2.5  ±0.1  mm.  Each  ice  wafer  was  then  examined  under  polar¬ 
ized  light  to  ensure  that  it  was  a  single  crystal.  They  were  then  stored 
like  the  soil-deuterium  wafers  until  they  were  ready  for  use.  We  selected 
single  crystal  ice  as  opposed  to  polycrystalline  ice  to  guard  against  any 
signal  contribution  from  liquid  films  between  ice  grain  boundaries.  Nu¬ 
clear  Magnetic  Resonance  (NMR)  readings  were  taken  on  samples  of  the  single 
crystal  ice.  No  NMR  signal  above  the  background  reading  was  detected  at 
-0. 5°C. 


Plonui  of  Symmetry  J 

for  Diffusion  )  ; 


A  test  tube  150  mm  long  and  19  mm  o.d.  Deuterium  Diffusion  Espenmen. 

with  a  rubber  stopper  was  selected  to  contain 
and  seal  the  samples  for  NMR  analysis.  A  16- 
-o.d.  by  5-mm-thick  plastic  disk  was  placed  in 
the  bottom  of  the  test  tube  to  serve  as  a  flat 

support  for  the  samples.  Next,  a  frozen  soil-  y777777g£^-\ 

deuterium  oxide  wafer  was  added,  followed  by  an  IwmjwL 

ice  wafer.  This  sandwiching  of  ice  wafers  be-  pi0n«»  of  SymmuryJ  ~ 

for  Diffusion  \  .  tjjjjjjjjji 

tween  soil-deuterium  oxide  wafers  was  repeated  .  . 

_  f  - 

until  the  complete  detection  region  of  the  NMR  -  - —  - — -  - 

analyzer  was  filled  by  16  wafers  (Fig.  1).  Two  77///////; 

of  the  layered  samples  were  made  (called  sample  ///////// 

A  and  sample  B  here),  and  one  sample  of  the  Mor-  ^  y 

in  clay  mixed  with  7X  distilled  water  without 
any  ice  layers  was  also  fashioned.  The  samples 

were  prepared  in  the  -20°C  coldroom.  Figure  1.  Schematic  o 

layered  diffusion  expe 

All  three  test  tubes  were  brought  out  of  iment  showing  ice  and 

the  coldroom  and  an  NMR  reading  taken  on  each.  8°^  wa^®rs  ®n<* 

of  symmetry  for  diffu- 

The  method  of  making  the  NMR  readings  and  sion. 

calculating  unfrozen  water  contents  from  the 

readings  is  discussed  by  Tice  et  al.  (1982).  After  the  MR  readings  were 
taken,  each  test  tube  was  placed  in  a  constant  temperature  bath  held  at 


Figure  1.  Schematic  of 
layered  diffusion  exper¬ 
iment  showing  ice  and 
soil  wafers  and  planes 
of  symmetry  for  diffu¬ 
sion. 


Table  1.  (MR  readings  on  layered  samples  and  blank  sample.* 


Time 

(hr) 

Sample 

A 

Sample 

B 

Blank 

Temperature 

(°C> 

■KWH 

29 

29 

51 

— 

0.28 

33 

32 

71 

-0.59 

0.43 

34 

34 

79 

-0.51 

0.83 

32 

35 

83 

-0.43 

1.33 

34 

35 

86 

-0.43 

1.83 

34 

35 

89 

-0.43 

2.83 

34 

33 

84 

-0.48 

18.3 

36 

36 

87 

-0.54 

22.3 

38 

37 

88 

-0.54 

26.8 

35 

35 

87 

-0.54 

43.0 

39 

38 

85 

-0.54 

50 

41 

41 

84 

115 

40 

85 

-0.56 

120 

41 

40 

87 

-0.54 

139 

43 

42 

87 

-0.51 

163 

44 

43 

87 

-0.54 

283 

45 

45 

87 

-0.54 

308 

46 

47 

88 

-0.54 

336 

47 

47 

87 

-0.54 

380 

47 

47 

88 

-0.54 

452 

47 

47 

87 

-0.54 

*Background  reading  was  26. 

Table  2. 
final  ice 

Deuterium  content  of 
and  soil  samples. 

Sample 

Atom  Z  Deuterium 

A-Ice 

5.54 

A-Ice 

5.44 

B-Ice 

4.77 

B-Ice 

4.75 

Avg.  5.13 

A-Soil 

30.24 

A-Soil 

31.77 

B-Soil 

27.01 

B-Soil 

32.55 

Avg.  30. 4 

-0. 5°C.  The  test  tubes  were  periodically  removed  from  the  bath  and  NMR 
readings  taken  over  19  days.  The  readings  are  given  in  Table  1. 

After  19  days,  the  sample  tubes  were  taken  back  into  the  -20#C  cold 
room  and  the  layers  of  frozen  soil  and  ice  from  the  two  layered  samples 


were  separated  and  put  into  discrete  sealed  containers*  Later,  both  the 
soil  and  ice  layers  were  analyzed  for  the  ratio  of  protlun  to  deuterium 
they  contained  using  an  Isotope  ratio  mass  spectrometer.  The  details  of 
this  analysis  are  given  by  Oliphant  et  al.  (1982).  The  results  are  shown 
in  Table  2. 


MATHEMATICAL  MODEL  OF  DIFFUSION  PROCESS 

At  the  equilibrium  temperature  of  -0.5°C,  all  the  water  in  the  clay 
disks  would  be  in  the  unfrozen  or  liquid  state  (Nakano  et  al.  1983).  As 
already  discussed,  the  diffusion  rate  is  about  5  orders  of  magnitude  higher 
in  the  liquid  state  than  in  ice.  We  therefore  made  the  assumption  that  in 
the  clay  the  ratio  of  hydrogen  to  deuterium  was  a  function  of  time  but  not 
of  position  in  the  clay  disks.  In  the  ice  disks  the  deuterium  concentra¬ 
tion  is  a  function  of  both  time  and  the  distance  into  the  disk.  In  the 
stack  of  alternating  ice  and  soil  disks,  each  ice  disk  is  sandwiched  be¬ 
tween  two  soil  disks.  This  gives  rise  to  a  plane  of  symmetry  at  the  center 
of  each  ice  disk  and  each  soil  disk  for  the  diffusion  problem.  To  model 
the  whole  diffusion  process,  only  one  half  of  a  soil  disk  and  one  half  of 
an  ice  disk  and  the  Interface  between  them  need  to  be  considered,  as  shown 
in  Figure  1. 

Letting  C  be  the  mole  fraction  of  deuterium  in  the  ice  phase  and  D  be 

2 

the  diffusion  coefficient  in  cm  / s,  we  can  write 

a2C  1  3C  n  ... 

17  '  Dir *°  (l> 

for  the  mole  fraction  of  deuterium  in  the  ice  phase  as  a  function  of  dis¬ 
tance  into  the  ice  x  and  time  t.  This  equation  holds  for  all  values  of  x 
between  the  interface,  x  *  0,  and  the  midplane  of  the  ice  disk,  x  *  l,  and 
for  all  times  t  >  0.  Because  of  symmetry,  we  can  also  write  at  x  ■  t  that 

£-o  .  <» 

A  mass  balance  for  deuterium  at  the  interface,  x  ■  0,  gives  the  equation 

30  t  JLl  2. 

3t  *  n2  i  l3xJx*0 

where  is  the  mole  fraction  of  deuterium  in  the  soil  disk,  nj  and  n2 
are  the  total  moles  of  hydrogen  plus  deuterium  per  cm  in  the  ice  and  soil 


phases,  respectively,  and  the  subscript  x  ■  0  means  that  this  derivative 
should  be  taken  at  the  interface.  We  can  also  write  at  the  interface 

<4> 

where  H  has  units  of  moles/cm  s  and  represents  the  rate  of  mass  transfer 
across  the  ice  disk/soil  disk  Interface  attributable  to  a  unit  concentra¬ 
tion  gradient.  The  higher  the  value  for  H,  the  lower  the  resistance  to  in¬ 
terphase  mass  transport.  A  lower  bound  for  H  can  be  estimated  as  follows: 
Assume  that  transport  across  the  interface  oust  be  by  diffusion  through  an 
air  gap  between  disks  with  maximum  thickness  of  0.01  cm.  The  diffusion  co¬ 
efficient  for  water  vapor  diffusion  through  air  is  on  the  order  of  0.1 
cm Is.  Using  the  ideal  gas  equation  and  the  vapor  pressure  of  water  at 

0°C,  we  can  calculate  that  the  maximum  concentration  of  deuterium  in  the 

7  3 

air  would  be  about  4x10“  moles/cm  .  We  can  calculate  H  to  be 

H  “  0.1  cm2/s  TTTrr -  x  4xl0“7  moles/cm3  =  4xl0“6  moles/cm2  s.  (5) 

u l  cm 


Smaller  values  for  the  air  gap  would  give  a  higher  value  for  H. 

Equations  1  through  4,  with  the  given  boundary  conditions  and  the  ini¬ 
tial  conditions  Cj i  “  C^0  and  C  =  0,  can  be  solved  using  Laplace  trans¬ 
forms  (Carslaw  and  Jaeger  1950).  The  final  solution  is 


and 


C 


a 2  -  KL  -a§T 
p  cos a  e  cos a8(l-x/i) 
s  s 


Cl 


C I 
1+K 


+  2KL  ‘ 


T 


where 

L  -  *h 
h  =«  H/niD 
K  *  nj/n2 
T  -  Dt/t2 


(6) 


(7) 


and  the  a8  are  the  positive  roots  of  the  equation 


and 


Ps  “  “s  +  (L  +  L  ”  2KL)Og  +  KL^( 1+K)«  (9) 

A  computer  program  was  written  that  calculates  C  and  C&  according  to  eq 
7  and  8  for  input  values  of  H,  nj,  n2,  D,  £,  Cj°  and  t.  This  is  given  in 
Appendix  A. 

RESULTS  AND  DISCUSSION 

Although  essentially  pure  D20  was  initially  added  to  the  Morin  clay 
samples,  the  samples  originally  contained  a  small  amount  of  water  and  also 
protons  associated  with  the  clay  matrix  that  would  equilibrate  with  the  ad¬ 
ded  deuterium.  The  initial  concentration  of  hydrogen  in  the  soil  sample 
can  be  estimated  as  follows  from  the  NMR  data. 

The  blank  soil  column  gave  an  NMR  reading  of  about  87  at  the  equilib¬ 
rium  temperature.  This  column  was  made  of  clay  mixed  with  distilled  water 
that  is  essentially  pure  H20  and  contains  only  the  natural  abundance  level 
(0.015%)  of  deuterium.  Subtracting  the  background  reading  of  26  from  87 
gives  61  counts  ascribable  to  the  hydrogen  in  the  added  water.  This  number 
must  be  divided  by  2  to  compare  with  the  readings  obtained  on  the  A  and  B 
samples  because  half  of ‘ the  volume  of  these  samples  is  made  up  of  single 
crystal  ice,  which  does  not  contribute  to  the  NMR  signal.  Thus,  a  signal 
of  30.5  plus  the  background  of  26,  or  56.5,  would  correspond  to  100%  hydro¬ 
gen  in  the  soil  disks.  The  final  NMR  readings  of  about  47  on  samples  A  and 
B  correspond  to  68. 8%  hydrogen  in  the  samples.  This  cotqpares  well  with  the 
mass  spectrometer  value  of  69. 6%  hydrogen  in  the  final  soil  samples  as 
shown  in  Table  2  (atom  %  hydrogen  *  100  -  30.4  »  69.6%).  All  the  1MR  read¬ 
ings  on  samples  A  and  B  can  be  converted  to  mole  X  hydrogen  values  by  simi¬ 
lar  calculations.  These  values  are  shown  in  Figure  2.  The  initial  values 
gave  NMR  readings  of  about  33.  This  corresponds  to  an  initial  hydrogen 
content  of  the  soil  disks  of  23X. 

The  ratio  ni/n2  used  in  eq  6  and  7  can  be  calculated  from  the  NMR  and 
mass  spectrometer  data  by  making  hydrogen  or  deuterium  material  balances. 
The  final  amount  of  deuterium  in  the  system  must  be  equal  to  the  initial 
amount  in  the  soil  disks.  This  leads  to  the  equation 

0.0513  ni  +  0.304  n2  -  0.77  n2  (10) 


> 


V  "J 
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Figure  2.  NMR  data  of  mole  Z  hydrogen  in  Che 
soil  wafers  as  a  function  of  time  (dashed  line 
represents  D  -  2.5x10”  0  cm/s;  solid  line  rep¬ 
resents  D  *  lxlO”10  cm/s). 

where  the  two  terms  on  the  left  represent  the  deuterium  in  the  final  soil 
and  ice  phases  and  the  term  on  the  right  represents  the  deuterium  in  the 
initial  soil  phase.  This  equation  gives  a  ratio  of  nj/n2  of  9.08.  Consid¬ 
ering  the  density  of  ice,  we  can  calculate  n^  to  have  a  value  of  0.102 
mole 8 /cm  ,  and  thus  n2  has  a  value  of  0.011  moles/cm  . 

Equations  6  and  7  were  solved  for  values  of  H,  n;,  n^,  i  and  Cjj,0 
equal  to  4xl0”6,  0.102,  0.011,  0.125,  and  0.77,  respectively,  and  for  two 

ig  2  in  2 

values  of  D,  1x10“  cm  /s  and  2.5x10”  cm  /s.  The  two  solutions  for  the 
mole  Z  hydrogen  in  the  soil  phase  are  shown  in  Figure  2  as  a  function  of 
time.  The  solution  with  D  *  1x10" 10  cm2/s  fits  the  data  well  at  times  less 
than  200  hours  and  the  two  solutions  bracket  the  data  at  times  from  300  to 
500  hours. 

Me  show  in  Figure  3  the  calculated  mole  fraction  of  hydrogen  as  a 
function  of  position  in  the  soil  and  ice  phases  at  several  different  times 
for  D  ■  lxlO"10.  It  can  be  seen  that  the  change  in  concentration  across 
the  interface  is  negligible  even  though  a  lower  bound  for  the  value  of  H 
was  used  in  the  calculation.  Thus,  interfacial  resistance  to  diffusion  ap¬ 
pears  to  be  negligible  for  this  system.  It  can  also  be  seen  that  times  on 
the  order  of  1000  hours  or  more  are  necessary  for  significant  approaches  to 
equilibrium  in  the  ice  disks. 
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Figure  3.  Calculated  concentration  of  hydrogen 
in  the  ice  and  soil  wafers  as  a  function  of  po¬ 
sition  at  several  times. 
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Figure  4.  Calculated  concentration  of  hydrogen 
in  water  and  ice  phases  for  a  hypothetical  exper¬ 
iment  using  ice  crystals  0.0025  cm  thick. 

In  Figure  4  is  shown  the  mole  fraction  of  hydrogen  as  a  function  of 

3 

position  at  several  times  for  ice  pieces  2.5x10“  cm  thick,  which  is  prob¬ 
ably  closer  to  the  size  of  ice  crystals  that  form  in  rapidly  frozen  soils. 
For  this  calculation,  ni  and  n2  were  set  equal  to  0.102  moles/cm  and  l  was 
set  equal  to  0.00125  cm.  All  other  parameters  were  the  same  as  for  the 
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calculation  shown  In  Figure  3.  The  interfacial  resistance  is  still  negli¬ 
gible  compared  to  the  resistance  to  diffusion  in  the  ice  phase.  Most  im¬ 
portantly,  at  this  size  for  the  ice  crystals,  equilibrium  is  approached  af¬ 
ter  only  a  few  hours. 

CONCLUSIONS 

Ue  found  that  using  a  value  for  the  diffusivity  D  of  deuterium  in  ice 
of  1x10“  cm  /s  in  an  analytical  model  of  an  experimental  diffusion  proc¬ 
ess  gave  a  good  fit  to  both  NMR  and  mass  spectrometer  data  obtained  for  the 
experiment.  We  also  showed  that  for  ice  having  a  thickness  of  0.25  cm,  a 
time  of  several  thousand  hours  was  required  for  equilibration  of  deuterium 
and  hydrogen  isotopes  between  the  ice  phase  and  surrounding  liquid  water. 
For  ice  thicknesses  on  the  order  of  2.5x10“  cm,  only  a  couple  of  hours  is 
required. 

A  mathematical  model  was  given  for  the  isotope  equilibration  process 
that  will  allow  an  estimate  of  the  required  equilibration  time  if  the 
thickness  of  ice  crystals  in  the  system  is  known.  This  can  be  used  to  help 
design  any  tracer  experiments  using  deuterium  in  frozen  soils. 
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APPENDIX  A.  COMPUTER  PROGRAM  THAT  CALCULATES  C  AND  C*. 


100  PRINT  ViNHU  HrNi  »N2»  »*l.i  rCO*- 

110  INPUT  H  f  N 1  ,  N2  r  1>  ,  l.  i  r  ('0 

115  PRINT  S41iH»Ni»N2rtt,Ll  ,00 

120  H1=H/N1/U 

130  L=L1*H1 

L  40  N=N1/N2 

ISO  FUZZ  10,1.0!- -.10 

160  DIM  A< 150) ,P< 150) 

l?0  FOR  1=1  10  150 

180  N=I 

190  Al=F‘I*<NH/2. 00000001) 

200  A2=PI*<N-l/2. 00000001 ) 

210  A3=<Al+A2>/2 

220  R 1 =T  AN  (  Ai  )  -I  *A  i  /  ( A 1  *  A  1  -K*L ) 

230  R2=TAN  ( A2 )  -L*A2/  ( A2*A2-K*L ) 

240  R3=TAN ( A3 )  -L*A3/ ( A3*A J-K*L ) 

250  IF  R3=0  THtN  310 

255  IF  A3=A2  THEN  310 

260  IF  SGN ( R2 ) "SON ( R3 )  THEN  2.90 

270  A1=A3 

2 SO  60  TO  210 

290  A2=A3 

300  GO  TO  210 

310  A(I)=A3 

315  PAGE 

320  PRINT  A3,R3 

330  NEXT  I 

340  DIM  X  <  51 ) , C2  <  51 ) 

350  FOR  J=1  10  20 
360  T=<J-1) *200000 
370  Tl  =  D*T/l_l/Ll 
330  31=0 

390  FOR  1=1  10  150 

400  P<  l  )=A<  I )'  4+(L*L+L-2*K*L)*ACO"2+K*!-*l  *<1+K) 

405  IF  A< I >"2*1 1>230  THIN  420 

410  31=S1  +  1/KU  >*F. *>'<-<  A (  T  )  '2 ) *  T 1 > 

420  NEXT  I 

430  C1.=C0/(l+K>+2*K*L*l..#00*6i 
440  PRINT  «?4i:T/3600,(Ji  ri-i'.l 
450  FOR  G-- 1  TO  6 
460  £2=0 

470  X < G  >  =  <  Q - 1 ) *L 1 /5 

480  FOR  1=1  10  150 

485  IF  A<  I )  "2*T  1>230  THI  N  510 

490  33=EXP  <  -  ( A  ( 1 )  '2 )  *  I  1  >  *603  ( A  <  I)  *  (  i  -X  <  (1  >  /L  l )  > 
500  S2=S2+(Aa>'"2-K*L)/<P()  )*00S<A<T)  )  >*S3 
510  NEXT  I 

520  C2<0)=C0/( i+K>-2*L*C0*S2 
530  PRINl  fe41JX<U>» 0201)  ,1-02(0) 

540  NEXT  0 
550  NEXT  J 


J'  <  V; 


f\'r 

YA 


i  w  - .  * . 


END 


FILMED 


1-85 


DTIC 


